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ABSTRACT 

A method to determine a bound of structural loads for a spacecraft 
mounted on a launch vehicle is developed. The method utilizes the interface 
shock spectra and the relative impedance of the spacecraft and launch 
vehicle. The method is developed for single^degree^of^freedom models and 
then generalized to multidegree-of»freedom models. 
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I. INTRODUCTION 

This report deals with the determination of structural design loads for 
a spacecraft mounted on a launch vehicle for which the range of significant 
resonant frequencies overlaps that of the spacecraft. Several methods have 
been used in the past to obtain spacecraft structural design loads. 

1 Analytically, the most elegant method to obtain loads has been to per- 

form transient analysis on a composite mathematical model of the space- 
craft and launch vehicle. This requires the availability of forcing functions, 
for example, transient chamber pressure of the vehicle engines, that are 
I applied to the composite model to compute accelerations and member loads 

for the entire composite structure. The main advantage of the transient 
[ I analysis method is that it leads to a lighter -weight design. However, the 

; I method is costly, leads to a complicated in. , face, and the results are sen- 

i I sitive to changes in the structures. Fabrication tolerances prevent exact 

\ definition of the model, and it is impractical to analyse all possible com- 

f \ binatlons of extreme tolerances. When this approach is used, testing must 

! I be limited to the derived loads. 

1 I Another widely used procedure has been to use the shock spectra of 

t 1 spacecraft/ launch vehicle Interface acceleration measured in previous flight 

I I define an environment more hostile than the actual tught. The envelope 

j of the shock spectra is used for the definition of an equivalent sinusoidal 

1 spacecraft base acceleration. Design and test have been made for either 

base control or response control. The advantage of this procedure is its 
simplicity, but the disadvantage is that it is sensitive to damping estimates, 
leads to heavy structural weight, and may not always be conservative. 

This report proposes a method for the determination of design loads 
that are realistic upper bound* of the flight loads. The method utilises 
shock spectra of launch vehicle accelerations in a rational manner and 
! 1 Introduces the relative impedance of the spacecraft and launch vehicle. It 

j uses global modal properties such as rigid elastic mass coupling which are 

less sensitive to minor modifications than the modal displacement of a pre- 
i 1 selected degree of freedom as for the traditional base acceleration approach. 
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The main objectives sought in the method are low cost of analysis, 
timeliness within the design schedule and high reliability, These objectives 
are attained at only a moderate expense of structural weight compared to a 
transient analysis design. A major advantage of the method is that the 
interaction between spacecraft and launch vehicle is introduced in a simple 
manner which permits a spacecraft loads analysis cycle to be performed 
without the necessity of a spacecraft/ launch vehicle composite modal analy- 
sis and flight simulation. Thus, the interface problems between the organi- 
zations responsible for the spacecraft and the launch vehicle are greatly 
minimized. Other advantages are low sensitivities to the details of the 
mathematical model of the structures, to damping estimates, and to minor 
structural changes. The method is recommended when rapid convergence to 
final loads is more important than the weight saving possible from detuning 
of spacecraft and launch vehicle modes. 

Finally, the method is a "follow on" method since it requires that 
spacecraft/l\unch vehicle interface flight accelerations have been either 
measured or calculated for previous flight using the same launch vehicle. 

n. METHOD 


A. SUMMARY 
1, Assumptions 

The assumptions on which the present method is based are; 

(1) The spacecraft and launch vehicle are adequately represented 
by linear mathematical models. 

(2) The shock spectra of the interface acceleration of previous 
flights are available from which an envelope can be constructed. 
The method is established for an unloaded Interface, but if no 
data is available for an unloaded interface, it is assumed that 
shock spectra for such an unloaded interface is within the 
envelope based on the loaded Interface data, 

(3) Response shock spectra of future flights will be within a con- 
servative envelope of those obtained from the prior flights. 
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(4) The frequency distribution of the source of the disturbance is 
relatively flat compared to that of the interface acceleration 
response. 

Note that most of the enumerated assumptions are also invoked for 
other methods of loads analysis. 

2. Procedure 

The method described here requires that the spacecraft to be analyzed 
and the launch vehicle be represented by their normal modes. Each space> 
craft cantilever mode and launch vehicle mode having resonant frequencies 
in the same range are paired. Then the largest possible spacecraft modal 
response is determined by allowing artificial shift of the two resonant fre- 
quencies to produce "tuning” between the two modes. 

Shock spectra of spacecraft/ launch vehicle interface acceleration of 
previous flights are needed to define the basic load. The relative impedance 
between launch vehicle and spacecraft modes is utilized to determine a 
reduction factor C of the dynamic response from that of a spacecraft with 
negligible impedance, i.e. , from the shock spectra S . 

The shock spectra is that of the unloaded interface (no spacecraft 
attached on the launch vehicle) and is assumed to have peaks, at the launch 
vehicle resonant frequencies, vdiose amplitudes are those of an envelope S 
of the available shock spectra. This envelope is assumed to enclose 

the shock spectra of the unloaded launch vehicle. 

The interface reaction and the member loads are determined for all 
possible reasonable pairings, and the corresponding absolute values of the 
loads are combined. The combination also includes the quasi- static case, if 
present. 
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Terminology 


The nomenclature used herein is presented as follows: 


(1) Symbols 
C 


c 



u 

Z 

V 

6 

e 

e 


reduction factor (maximum of 
defined in Eq. (21) 
damping 

reaction (positive on lavinch vehicle) 
equivalent applied force at interface 


spring stiffness 
mobility 
modal mass 

effective mass (modal mass normalized to be equal 
to elastic /rigid body coupling) 

shock spectra for damping 

time 

modal displacement 
displacement 
impedance (1/M) 

2 c/c^ (damping ratio) 

(w/w^ - w*'/w)/y*' 

(w/o)^ - wVw)/Y 





i 

T 

w 


modal motion of interface 
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u» circular frequency 

[3 a rectangular matrix 

I I a vector 

,2, Sub.cripts and Superecript. . Subscript, are u.ed to Indicate 

degrees ot freedom both physical and generalised. Superscripts 
are used as modifier, of the symbol if any ambiguity remain, 
after consideration of any subscripts. 

c critical 

D data point 

E equivalent 

e envelope 

I interface 

i mode of spacecraft 

k mode of launch vehicle 

L launch vehicle 

o initial value (without spacecraft) 

R relative or reduced 

S spacecraft 

1 8ingle-degree-of-freedom oscillator interface 


(3) Operators 


derivative of ( ) with respect to time 

equal 

approximately equal 


of terms for each value of i 


= definition 

X 1 absolute value of x 
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tranf^pose of C 1 
inverse of C 3 


B. SHOCK SPECTRA 

1 . Definition 

The shock spectrum S^(w) is the maximum acceleration response of 
a single-degree-of-freedom oscillator of mass m, damping y - Zcfc^, and 
natural frefiuency w to abase acceleration u^ (T) (Fig. 1). 

SV) » |V'^' 

where is the acceleration of the mass m relative to the base. Here 

the base is the spacecraft/ launch vehicle interface. 

It is noted that: 

(1) The base acceleration *d°(T) is an input unaffected by the pres- 
ence of the oscillator; i. e. , the impedance of the oscillator is 
infinitesimal compared to that of the base. 



Fig. 1. Shock spectrum 
oscillator 
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(2) The shock spectrum is a function of two variables, the natural 
frequency u> and the damping y of the oscillator, and exhibits 
a series of peaks of various heights and sharpnesses. 

(3) The maximum reaction force on the spring connecting the 
mass m to the base is equal to the inertial force of the mass m. 

F j = m s''’(w) (2) 

Therefore, the shock spectrum is a measure of the reaction force on 
a unit mass. This observation is an essential part of the method presented 

in this report. 

2. Envelope 

Because in structural design one has to consider the most critical 
loading condition, it is customary to envelope the ensemble of shock spectra 

[s^m]". n=l,2 N. of the interface acceleration [u^(T)J of previous 

flights, measured or calculated, using the same launch vehicle in order to 
determine a shock spectra envelope which is substituted for S^(w). The 

construction of the envelope introduces a degree of conservatism that 
depends upon how it is done. Figure 2 shows the traditional enveloping tech- 
nique by which straight lines are drawn tangent to the highest peaks of the 



log ui 

Fig. 2. Shock spectra 
and envelope 
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ensemble of all shock spectra or higher than those peaks depending upon the 
desired degree of conservatism. This envelope has the advantage of being 
simple to construct. However, it is very conservative since it has the dis- 
advantage of flattening out all peaks of the shock spectra and therefore 
ignoring valuable chai acteristics of the source of the base acceleration 
u°(T), namely, the peaking effect in terms of frequency. This effect will 
be reintroduced in the method described. 

C. USE OF THE SHOCK SPECTRA 

The shock spectra envelope S^® can be utilized to determine the 
spacecraft/ launch vehicle interface load and also the loads in the entire 
spacecraft. The oscillator (Fig. 1) can represent each cantilevered mode of 
the spacecraft. Calling m" the effective mass of a mode and y" the damping 
of that mode, then the maximum reaction force on the base of the spacecraft 
for that mode can be written as: 


where S^® is the value of the shock spectra envelope corresponding to the 
damping\^ and the natural frequency of that mode. Note that Eq. (3) 
relates the loads directly to the shock spectra and with the proper damping, 
and therefore represents a significant improvement over the traditional 
method of relating loads to static or sinusoidal accelerati on. 

However, Eq. (3) can be further improved by noting (1) that the enve- 
lope S'^® ignores the peaking character of each shock spectrum and (2) that 
the impedance of the spacecraft is not infinitesimal compared to that of the 
launch vehicle, The relative impedance between the spacecraft and launch 
vehicle reduces the interface motion and the reaction on the spacecraft and 
also shifts the frequency at which the maximum loading occurs. 

Therefore, a reduction factor C S 1 is introduced which depends on 
the shape of shock spectra peaks as well as the relative impedance, and 
Eq. (3) is modified as follows: 


F 


R 

I 


= S^® C 


(4) 
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The superscript R has been introduced as a reminder that this reaction has 
been reduced as a function of the relative impedance and peak shape. The 
object of the next section is to determine the reduction factor C. 

D. DETERMINATION OF THE REDUCTION FACTOR 


1. Equivalent Forcing Function 

To find the reduction factor C it is convenient to introduce a hypo- 
thetical sinusoidal forcing function F^e^^*^ applied at any convenient point 
E of the launch vehicle (Fig. 3), which for a given frequency u produces an 
interface motion with the same shock spectra as for the flight data. 

An expression of F^ proportional to is sought. Let us assume 
that F^ is applied to the unloaded launch vehicle, i.e. , without the space- 
craft (Fig. 4). 

Using the mobility concept (Ref. 1) for the unloaded launch vehicle, 
we find that the interface velocity | Ujj, | is the mobility | sit the inter- 

face due to a load at point E times the load F^ at point E. 




(5) 


The mobility | | of the launch vehicle is the sum of the mobility 

due to each of its modes k. In terms .^f modal frequency nodal 


SPACECKAFT / , \ 


LAUNCH VEHICLE 



Pf 


Fig. 3. Composite structure 



Fig. 4. Unloaded 
launch vehicle 


i 

i 
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mass m^, modal damping y * niodal displacements ♦£, and the 
frequency u» we have 




k k k 
w m Y 



(6j + 1) 


( 6 ) 


where 



Equation (6) exhibits peaks that occur at each resonant frequency of 
the unloaded launch vehicle. Assuming that the source of the disturbance 
has a relatively flat frequency distribution, we find that the shock spectra of 
the unloaded interface accelerations will also have peaks at the same 
resonant frequencies. The envelope of all shock spectra will be con- 

trolled by the amplitude of those peaks. This observation is fundamental for 
VC 

the handling of S 

Let us now note that in practice, flight measurements are most likely 
made only foi a loaded interface, i.e. , for the launch vehicle carrying a 
payload. However, since the payload has a small mass compared to the 
launch vehicle, the presence of this payload does not appreciably change the 
amplitude of the peaks of the unloaded launch vehicle interface acceleration 
shock spectra (Appendix A), although the frequencies of those peaks are 
somewhat shified. In addition, secondary peaks will occur near the canti- 
levered spacecraft frequencies, but these peaks will be lower than those due 
to the launch vehicle and hence will not control the shock spectra envelope 

It is therefore postulated that the shock spectra envelope Is con- 
trolled only by the amplitude of the launch vehicle peaks and is unaffected 
by the presence of the spacecraft; i.e. , the shock spectra can be 

determined from flights with a loaded interface. Note that the postulate is 
not valid if the data comes from a flight where a launch vehicle resonant 
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fr6qu6ncy is tuned to & spececraft resonant frequency. In this case, the 
controlling launch vehicle peak will have a lower amplitude than for an 
untuned case (Appendix A). However, this tuned case is unlikely to occur 
for all spacecraft modes of all prior flights; the envelope will be determined 
by the untuned case or by other flights without tuning. 

It was noted earlier in Section II-B-2 that the envelope S^® flattens 
out the peaks of the real shock spectra [s^(w)]”. In order to reintroduce the 
peaking characteristic of the real shock spectra, a h ypothetic al curve S(Y»®) 
enveloping any is defined as S(Y,o) = + l for 6 or i 1 

(Fig. 5). This enveloping curve has the amplitude of the peak at the level 
I S"^® of the simple envelope but allows the peak to shift in frequency along 

\ s''^®. The coefficient a controls the sharpness of the enveloping curve and 

f is chosen such that all peaks of each available shock spectra data are 

I enveloped. The limitation of 6^or < 1 corresponds to the half-power point 

■ cutoff to account for close peaks and background noise that would otherwise 

create problems in implementation. 

The amplitude Fj. of the equivalent forcing functions will now be 

determined. To relate F-, to the flight acceleration shock spectra envelope 

" — i 

we attach to the unloaded interface an infinitesimal oscillator of mass m , 

natural frequency w', and damping Y^ (Fig. 6). 

The maximum amplitude of the velocity response Uj of the mass of the 
i oscillator due to F^ occurs at and is given by 

! 1, = • wJ-j-Fg/jY* (7 



Fig. 5. Hypothetical curve enveloping peaks of shock spectra 
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Fig. 6, Infinitesimal oscillator 
on unloaded interface 


The amplitude of the equivalent force Fj, is obtained by equating u^ with 
the shock spectra envelope modified by the assumed shape of a shock spectra 
peak. We have 


u. s S 

1 




/j 




+ 1 


( 8 ) 


Since this envelope depends only on the amplitude of the peaks, one 
can consider an interface motion due to only one launch vehicle mode at a 
time. Then combining Eqs. (7) and (8) the equivalent force Fj, is 

p ( 9 » 

where is the mobility corresponding to one launch vehicle mode in 
Eq. (6), in which the arbitrary location of point E has been taken at the 
interface I for convenience (#£ ■ The equivalent force Fj. is an 

external force at the interface I , 




I 


m y w (8j + 1) 


( 10 ) 
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2, Interface Reactton 

When the launch vehicle is loaded by the spacecraft (Fig. 7), a reaction 
force |Fj| is induced by the epaceoraft on the launch vehicle interface. 

For the launch vehicle the interface velocity | | i® mobility 

{*^IE } ** interface I due to a force at E times the force F^, at E, plus 
the mobility [mJj ] at the interface due to a force at the interface, times 
the reaction force | Fj | at the interface I; 


{ Uj I = 



( 11 ) 


For the spacecraft the interface velocity is the mobility [ Mjj J at the 
interface due to the force at the interface, times the reaction | -F j J on the 
spacecraft: 


lii(= [Mfj]{.r,} 

Combining Eqs. (11) and (12): 


(12) 


(13) 


SPACfCKAFT 


A 


LAUNCH VIHICLI 



•f,, 6, 


Figi 7. l-sunrh vehicle/ spacecraft 
composite structure 
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Rearranging Eq. (13) and using the superscript R as a reminder that 

+ (14) 

In Eq. (14) the mobilities are three-dimensional and contain all the 
modes of the spacecraft and the launch vehicle. However, the essence of 
the method leading to the reduction fadtor C can be presented at the single- 
degree-of-freedom level. Thus, the development of the method will utilize 
single-degree-of-freedom models, one for the spacecraft and one for the 
launch vehicle, with the data point at the interface and along the line of 
motion (Fig. 8). This special case will then be generalized. 

Ic i 

The modal mobilities M for the launch vehicle, Eq. (10), and M for 

L S 

the spacecraft will be substituted for and respectively, in Eq. (14). 

The modal mobility is the sum of the mobility of the effective mass 
and the mobility of the modal spring damper system connecting m^ to the 
interface (Fig. 9). 


M 


i 


1 

. — i 
j m u> 


+ 


1 

~rrj7~ 

c + k /wj 


(15) 



Fig. 6. Modal model of composite structure 
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Jj 
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Fig. 9. Relative velocities for 
cantilevered modes 


The effective mass f5l' is determined as the ratio of rigid/elastic 
mass coupling scjuared to the modal mass (Ref. 1). 


= JL |mjJ^ {"'ll} ^ 


( 16 ) 


m 


Expressing the damping and stiffness of Eq. (15) in terms of modal 
frequency damping ratio Y^, and effective mass m\ we obtain; 


= 


h> 




(17) 


Introducing e* (w/w^ - wVw)/y^ in Eq. ( 17) and neglecting y^w/wS 
as compared to unity we finally have 




( 1 &) 
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In addition, for single-degree-of-freedom models, Eq. (14) reduces to 




— 1 k 


(19) 


Substituting Eqs. (9), (10), and (18) into Eq. (19) and considering the 
absolute value, the following expression is derived for the reaction force: 


_R — i eV® 
Fj = m S' 




( 20 ) 


where 


i 



U) 



6 ^ + 1 


e + 1 - 6e)^ + ( 6 + e)^] (6^0 + 1) 


(21) 


with 



( 22 ) 


4 

w m 


^i^k 


For structural design, we are interested in the maximum value that 
the reaction can reach. This reaction depends upon two parameters, 

6 and E, for given vSv*', 0, and or. The two parameters 5 and e measure 
the proximity of the two resonant frequencies w of the spacecraft and 
the launch vehicle and the dummy frequency w. 
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A conunon approach for lighter -weight design is to avoid tuning between 
the two resonant frequencies and However . designing for well- 

separated resonant frequencies requires costly, accurate structural analysis 
and is difficult to achieve because of \inpredictable variation in fabrication 
tolerances. The method proposed here is to provide for the tuned condition 
even if tuning will not exist in practice. The method, although not optimum 
for weight, is cost-effective. It also increases reliability since only simpli- 
fied structural analysis is required and the design loads are for inadvertent 
tvining that may occur from resonant frequency shifts due to fabrication toler- 
ances. Heuristically, it is expected that the maximum of the reaction F j 
occurs for = w*' = w. However, close examination of Eq. (22) shows that 
the maximum of occurs for a condition somewhat different from that of 

tuning. The maximum value C of the coefficient 9^ in terms of the two 

i k 

independent variables 6 and e (correspondingly /w andw /w) has been 
determined numerically for convenience. Therefore, Eq. (20) becomes 




F^ < C 

C = C (Y^. e^a) (23) 

1 

e = — e 

The values of C as a function of v''. and 6^ are shown in Fig. 10 for 
values of a given in Fig. 11. 

3. Typical Values of a 

Typical values of the parameter a defining shock spectra peak shapes 
for the Titan IIIC booster are given in Fig. 11. The accelerations from 
which the shock spectra were obtained are at the top of the Centaur stub 
adapter. There is a great deal of scatter In the data, but a trend toward 
sharper peaks with lower damping is evident, A curve whose slope is 
approximately that of the trend is also shown. The analytic expression for 
that curve is 
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No attempt has been made to find a trend with respect to frequency. 
Values of a greater than 1/6^ will not be utilized as indicated in 
Section II- D- 1 , 

4. Typical Values of C 

The approximation to a is used to evaluate values of C which can be 
used for payloads supported on the Titan IIIC booster. The values of C 
shown in Fig. 10 indicate that C is more sensitive to the geometric mean 
Y = V of spacecraft and launch vehicle damping than to either one 

separately. The reduction effect is greater for light damping, but the shock 
spectra is expected to be higher for light damping. Thus, the two effects of 
damping oppose each other, and the reaction is not overly sensitive to 

damping. 

5. Ty pical Shock Spectra Envelopes 

Shock spectra envelopes for axial and lateral acceleration at the top of 
the ?cub adaptor of the Titan IIIC are shown in Fig. 12. These envelopes are 
based on simulated flights using engine chamber pressures to generate forc- 
ing functions and both Viking and Helios spacecraft models. 

6. Damping 

If the spacecraft damping is not known accurately, the reaction can 
still be bounded by evaluating the reaction using both a low and a high esti- 
mate of the damping in combination with the corresponding shock spectra. 

E. GENERALIZATION 

1. Modal Reaction Force 

The derivation of the reaction force F^ in Section II-D-2 was done for 
a single-degree-of-freedom oscillator along its line of motion. The 
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Fig. 12. Shock spectra envelopes at 
top of Titan III-C 

derivation will now be generalized to a rnultiaimensional interface I to deter- 
mine the multidimensional reaction force j F j | at this interface. 

Call ] ffij. I the mass coupling between a spacecraft cantilever mode of 
effective mass m. and the unit rigid body interface motions; then the equiva- 
lent interface motion u^. (Fig. 8) along a single-degree-of-freedom repre- 

1 I T „ . 


sentation of the mode i is(l/mi)|Fnj 
I Uj I at the interface I. 


(Ref. 2), times the real motion 


i I 1_, I T j I <241 

ri| 

The reaction vector |f|| at the I interface is ( 1 /mM{rfTj^{ times F^ pre- 
viously calculated, along the line of motion of the single-degree-of-freedom 


system 
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The motion juj | at the interface I is the mode shape |of that interface 
for mode k times the modal motion of the launch vehicle 

k 

The motion at the data point is the mode shape evaluated at the data 
point times the modal motion of the launch vehicle U^: 


k . k ..k 

“d = *D W 


From Eqs. (26) and (27) the motion at the interface is 




T j shock spectra peaks and envelopes can be scaled by the same 

factors since they are due to launch vehicle modes, 




This allows us to use flight data that has been recorded at points other than 
the interface. Analogous to Eq, (24), 




Combining Eqs. (23), (25), (29), and (30), 


I ' ■jj? {"li } "‘'Sf {"'M I 
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which simplifies to 




m 


( 32 ) 


where the bar has been omitted for m. and |mj. | since the triplet 

is independent of modal normalization. The normalization for an effective 
mass m., if used, should be applied to both the mass coupling and the 

mode shapes U'lof the spacecraft. The square of the factor is applied to 
m' so its effect on the reaction and later on the displacement vector cancels. 
This is expected since the method of normalization cannot affect the results 
but is a useful aid in interpreting the model. 


For parallelism with Eq. (23) the follo\^dng definitions are made: 


H] * 


I T 

lil \^li\ 


(33) 


lc>^ke I 

I I I 


■if-KI 


(34) 


Substituting Eqs. (33) and (34) into Eq. (32), we obtain 


{Ff} = } 


(35) 


iri*'! are a 


Note that Wjj and {*^1 ) 
to the number of degrees of freedom at the interface 


matrix and a vector with order equal 
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Z, Evaluation of the Reduction Factor C 

The reduced reaction has the same form for the single-degree-of 

A 4 U 

freedom models of Section Il-D-2. To evaluate C , curves of Fig. 10 are 
used for which the coefficient 6^ becomes 6^^ as follows. From Eqs. (24) 
and (26) we have: 


m ' ' 


(36) 


ic 1 k 

The part preceding U in Eq. (36) is the mode shape that 

replaces of Section II-D-2. 


* 




(37) 


Substituting Eq. (37) in the definition of 6 , we obtain 


e 


ik _ 

SB 



(38) 


which simplifies to 


Ik j ,59) 

w m Tti y y 

Note again that the renormalieation of the modes to obtain the effective 
mass ffij is not required in Eq. (39) as was the case for Eq. (37). 

3. Modal Mem ber Loads 

Member loads | pM due to each mode are obtained by premultiplying the 
displacement vector juM of the spacecraft by the stress coefficient matrix 
CSC 3 of the member. 
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(40) 


1p‘1 ” Csc3)o‘l 

The displacement vector {u'| is obtained by multiplying tiie spacecraft 
modejT^i normalized for equivalent effective mass by the relative motion 

- u*^ - u (Fig. 9) of the effective mass representation of the spacecraft 
R i I 
mode. 

A conservative estimate of member loads for spacecraft mode ^^d 
launch vehicle mode k is obtained by replacing Uj^ by ii. = Sj^ - Uj . 

The shock spectra at the data point is che sum of the data point acceleration 
h* and the relative acceleration Uj^ of the oscillator. The contribution of 
each is not known since such information is lost in the total acceleration 
shock spectra processing. However, since the interface base acceleration 
will be used without reduction while the modal acceleration will be multiplied 
by a factor less than 1.0, overestimates of the base acceleration contribution 
are conservative (see Appendix B). The maximum value of the data point 
acceleration as measured by the tail of the shock spectra at frequencies 
greater than those of interest will be used as the estimate of the interface 
base acceleration ti. in the same direction. It is also conservative to 
assume that the rigid body and modal acceleration are in phase. Therefore, 
it is conservative to estimate the relative acceleration as the data point 
shock spectra less any interface base acceleration of significance. The 
relative motion vector i uM is then estimated as: 


The total motion uj'* in Eq. (41) is obtained from the acceleration 
estimate uj'* of the effective mass of the spacecraft and the reaction 

force 


k* 

i 



w 


S 



T 

U> 


(42) 
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Using Eqs. (25), (32), and (42), we have, at w = w , 


^ik /^Yc . y max\ 

z\S-lh—l irri \T/ 4 ,k) 

"i ' /\\2 7 \l) 


( i)^ .k _i 
\w / ♦_ Tf\ 


D 


(43) 


Combining Eqs. (40), (41) and (43) we obtain the member loads 

-ik /eve max\ 

{P*}- 


(44) 


where again the bar has been omitted. 

4. Total Loads 

In the previous sections only the modal loads for one pair of modes are 
calculated at a time. Consideration will now be given to the combinations of 
the modal loads to obtain a bound of the total loads. This combination pre- 
sents a difficulty for the present method. When excited, all the modes are 
not responding "in phase," so load maxima are not occurring simultrineously 
for ah modes. However, since the modal loads developed in the previous 
sections are bounds, the phase Information is not available. 

A bound of the total load can be readily calculated by assuming that at 
some point in the time history all the modal loads in a given direction will 
add. That is. 


■'f ^ i: kr 


in one direction 


(45) 


and 


P < 


£|,. 


one direction 


(46) 
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This bound is overly conservative unless only one pair of modes is predomi- 
nant over all others. 

If two or more pairs give similar modal loads, their relative occur- 
rence should be examine*!. This would be particularly true If the frequencies 
of the pairings are widely separated. Since the excitation is transient in 
nature, one maxima is most likely to occur well before the other. 

The conservatism built into Eqs. (45) ard (46) can be greatly reduced 
by introducing some elementary statistics in the formulation for the combina- 
tion of modal loads. In Eqs. (45) and (46) it is assumed that all modes of the 
spacecraft are paired and tuned simultaneously with a mode of the launch 
vehicle. This situation, although theoretically possible, is actually quite 
unlikely in practice since only a few, if any, spacecraft modes will be tuned 
to the launch vehicle modes. Consequently, one can treat all the bounds of 
the modal loads as statistical quantities and introduce the root sum square 
(rss) technique to evaluate the loads as follows: 


and 




EC'-")' 


with high reliability 



wiUi high reliability 


(47) 


(48) 


An alternative method to subtracting base motion from total acceleration 
shock spectra would be to generate relative acceleration shock spectra 
directly and include rigid body accelerations end the associated quasistatic 
stresses in the RSS summation. This has not been done, since only a total 
motion shock spectra evaluation program was on hand. 

In the problem used as an example, a spacecraft supported on a launch 
vehicle, ttie only significant static acceleration is along the vehicle axis, and 
the data point in that direction is either at the interface or has almost the 
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same modal displacement in all modes of interest. For such a problem the 
method described has the advantage over the alternate method that greater 
assurance of bounding is possible with negligible penalty. 

Let us now make some additional remarks concerning each of the 
modal loads. Equations (35) and (40) give the desired interface reaction and 
member loads for mode i of the spacecraft when the motion at the interface 
is predominantly due to mode k of the launch vehicle. Use of the launch 
vehicle mode which is the most hostile is conservative. If more than one 
is suitable to use, the one resulting in the lowest reaction should be 
used, as all are bounds. 

The lumped mass model of each normal mode represents the rigid 

body inertial properties as well as the elastic properties. Quasistatic accel> 

eration is represented even for cutoff transients if the quasistatic boundary 

conditions are included in the derivation of the shock spectra. However, 

only the dynamic portion of the response is subject to reduction by the fac> 
ilc 

tor C . The appropriate quasistatic acceleration should be subtracted 
from the shock spectra before it is used. Also the effect of the quasistatic 
load must be added unreduced to the summation of the contribution of each 
mode. 
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APPENDIX A 


VARIATION OF INTERFACE ACCELERATION IN TERMS 
OF SPACECRAFT/LAUNCH VEHICLE INTERACTION 

The effect of prior spacecraft on vehicles from which acceleration 
data was taken is now inspected. Because of the assumption made in Sec- 
tion II-D-1, the amplitude of the effective forcing function at the launch 
vehicle resonance without a spacecraft is the same as at the composite reso- 
nance of a launch vehicle with a spacecraft attached. The ratio of maxi- 
mum motions at the interface unloaded and loaded is the ratio of unloaded to 
loaded maximum mobilities 


= 




(A-1) 


The impedance of the loaded vehicle is the sum of the impedance of 
the spacecraft and the launch vehicle 


M 


ik 

11 


M 


11 


M 


11 


(A-2) 


The mobility of the unloaded vehicle is (see Section II-D-1) 




k k „k , 

w m V ( 6.1 + 1 ) 


(A- 3) 


The unloaded launch vehicle has a maximum mobility when it# imagi- 
nary part is zero; i.e. , 6 - 0 . 





k k^k 
u» m Y 


( A-4) 


j 
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The mobility of the spacecraft is (sec Section II-D-l) 


M 




(A-5) 


w m 


Therefore. Eq. (A-2) becomes, after reduction 

[{1 + e^)6- 0e] j 


M 


u (M^r 


1 + 


e 

+ 1 


+ 1 


(A-6) 


The loaded launch vehicle mobility is a minimum for a zero 

imaginary part in Eq. (A-6): 

(I + e^) 6 - ee = 0 


(A-7) 


Substituting Eq. (A-7) into Eq. (A-6) we have: 



(A-8) 


The terms e and 6 are not independent as both contain w. The term 
r « will be taken as the Independent variable. 

Rearranging Eq. (A-7). we have 



For » 1 or results insensitive to e, 

6c a 6 


(A-9) 


(A-10) 
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Substituting the definition of e , 6, and 6 into Eq. (A- 10), we have 





u) m I ^ j I 
^ k' i 

w m V Y 


Rearranging Eq. (A- 11) we obtain 


Let 


Then 



n^-(r + p + i)n+r = o 


Solving for Jl, we obtain 






i)^ - 


4r 


From the definition of e and n, 


c = 



(A-11) 


(A-12) 


(A-13) 


(A-14) 


(A-15) 
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Finally, from Eqs. (A-1) and (A- 8), 


Mi 


max 


e 




max 


e + 1 


(A-16) 


Equation (A-16) is the required expression of the ratio M . Fig- 

ure A-1 shows the ratio .5? as a function of r for Y /V =1 and several Y 

and P. Figure A- 2 shows the ratio oi x for y = 0* 02 and several p and 

y'/Y^. Figure A-3 shows the ratio M as a function of 6^ for r = 1 and 
i Ic 

several Y /Y and Y. 

k i • 

Those curves show that, except for the tuned case w = w , the ratio 
remains relatively constant, slightly higher than unity; i. e. , the amplitude of 
the unloaded interface acceleration is only slightly higher than that of the 
loaded interface. 
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Fig. A-2. Interface acceleration ratio, 










af^fendix b 


CONSERVATIVE REPLACEMENT FOR 
RELATIVE DISPLACEMENT 


It will be shown that, given C, 

0.0 < C < 1.0 


The following inequality holds: 

u* + C ( |«i 1 - - |“i ^ ^"i ' \ 

Where the base motion uj is replaced by an upper bound u* in phase 
with Uj , 


Two basic cases will be considered. The first provides for u|^. with 

the same sign as uj; the second provides for u^. with an opposite sign to 

u”. For convenience, uj will be assumed positive. This is done without 

iJss of generality by selection of an appropriate coordinate system. Let us 

call Un the reduced relative displacement 
R 

H C (U- - uj) 


Case 1: Uj^ > 0 

If u^ > 0 then > 0 

If u^ and u“ are both>0, then u^ is >0. 

The absolute value signs in the inequality to be proved can be removed: 
u* + C (u“ - u*) 2 Uj + C (u.® - uj) 
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Simplifying and rearranging, we obtain 


(1.0 - C) (u* - up > 0 

Since both factors are positive by definition, the inequality is proven in 
Case 1. 

CASE 2: < 0 

K u^ < 0, then u|j^ < 0 

Case 2 will be further subdivided into two subcases. 

SUBCASE A: uj > u^ 

K u° > u^, the absolute value on the right side of the inequality to be 
proved may be removed. 

u=:= + C { |u.° I - u*) i u“ + C (u.° = u[) 

If the inequality holds with the absolute value removed, which is 
Case 1, the inequality is proven in case 2A,, 

SUBCASE B: uj < u^ 

If u,° < u„ the absolute value sign on the right side of the inequality 
1 R 

to be proved may be removed after changing the sign of the terms included. 

u* + C (| u.° I - u*) > - uj - C (u° - up 

Rearranging, we obtain 

(uj + u*) (1 - C) + C (|u*| + up 2 0 

Both factors of both terms are positive by definition, so the inequality 
holds in Subcase 2B. 
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